Electron-Phonon Interactions in C28-derived Molecular Solids 
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We present ab initio density-functional calculations of molecular solids formed from C28-derived 
closed-shell fullerenes. Solid C28H4 is found to bind weakly and exhibits many of the electronic 
structure features of solid C60 with an enhanced electron-phonon interaction potential. We show 
that chemical doping of this structure is feasible, albeit more restrictive than its Ceo counterpart, 
with an estimated superconducting transition temperature exceeding those of the alkali-doped Ceo 
solids. 
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There continues to be intense research efforts in 
studying fullerenes for their diverse properties, not the 
least of which is their unusually high superconducting 
transition temperatures T c (up to 40 K in Cs3Ceo)i 
Theoretical studies show that many of the phenom- 
ena associated with this class of materials can be ex- 
plained within the electron-phonon mediated picture of 
superconductivity! 2 *!^^ The high T c of these materials 
relative to that of intercalated graphite is attributed to 
the curvature of Ceo' 3 i 7 i 8 Hence, solids based on smaller 
fullerenes possessing an electronic structure similar to 
that of alkali-doped Ceo may exhibit a T c enhancement. 

The goal of the present work is to identify new molec- 
ular solids analogous to the alkali-doped G§q solids 
with similar electronic structure and increased electron- 
phonon coupling. We consider solids composed of 
C28-derived closed-shell molecules (C24B4, C24N4, and 
C28H4) which can potentially have a large tunable 
density-of-states (DOS) arising from the narrow weakly 
broadened bands. Theoretical studies of the isolated 
molecules indicate that they should be stablej^iiSiiiiiS 
however, to our knowledge there are no studies of their 
solid forms. Here we present results of ab initio pseu- 
dopotential density-functional calculations to determine 
the structural and electronic properties of the solids in 
pristine and doped forms. There have been ab ini- 
tio studies of solids based on C28 (Refs. Il3lll4lll5l) and 
C :Ul ,i^J^i which exhibit increased electron-phonon cou- 
pling. However, these molecules form covalent solids 
differently from Ceo- We find that C24B4 and C24N4 
also form strongly-bonded solids, but that C28H4 exhibits 
many of the salient features of solid Ceo- Several doping 
scenarios are investigated for alkali doped C28H4 crystals, 
and a promising candidate is identified for high T c . 

First proposed by Krotoji^ the C28 molecule obeys the 
Td point group. This fullerene is produced with an abun- 
dance nearly as great as that of Ceo m the laser vapor- 
ization of graphite^ It is a very reactive molecule with 
dangling bonds localized on each of the four apex atoms 
shown in black for each C28 in Fig. 1(a). Because the 
C28 molecule is a tetravalent superatom, a reasonable 
candidate for a solid is a four-fold coordinated diamond- 



like structure. The C28 molecules can be four-fold co- 
ordinated with respect to the tetrahedron apexes as in 
Fig. 1(a) apex-bonded or with respect to the tetrahedron 
faces as in Fig. 1(b) face-bonded. Since the C28 molecules 
have dangling bonds localized on the tetrahedron apexes, 
the lower energy structure will be the former. Previous 
ab initio calculationsi^i^ have shown that the covalent 
bonds of the apex-bonded C28 hyperdiamond are very 
strong compared to the weak forces responsible for bond- 
ing in solid Ceo- 




FIG. 1: Schematic diagrams of the first nearest neighbors 
in the hyperdiamond structures (top) with examples of pairs 
of C2s-derived molecules shown below. The lattice sites in 
the solid are highlighted by the grey spheres enclosed in the 
tetrahedra. The two distinct orientations of the constituent 
molecules are represented by those of the tetrahedra. The 
nearest-neighbor pairs show two distinct bonding configu- 
rations: (a) apex-bonded C28 hyperdiamond forms covalent 
bonds (dotted line) between apex atoms colored in black; 
(b) face-bonded C28H4 hyperdiamond forms weak bonds be- 
tween six-membered rings with hydrogen atoms depicted as 
smaller black spheres. 
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In search of suitable molecular solids, we study hy- 
perdiamond structures constructed from three distinct 
closed-shell molecules - C24B4, C24N4, and C28H4, us- 
ing ab initio pseudopotential density functional meth- 
ods within the local density approximation (LDA). The 
SIESTA cod o 2 ?' 21 was used to perform conjugate gra- 
dient minimization of forces and stresses to find re- 
laxed structures. Though other possible solid structures 
for C28 have been proposed^ they are reported to be 
higher in energy than apex-bonded hyperdiamond and 
so were not considered in this study. Although all three 
molecules are stable, as was found previouslyj2*i2*ii*i£ we 
find that in the solid the C24B4 molecules are unstable 
and break apart. Of the two geometries investigated for 
C24N4, apex-bonded hyperdiamond is much lower in en- 
ergy forming a covalent solid like C28- The face-bonded 
C24N4 hyperdiamond is higher in energy and in addition 
shows significant hybridization in the conduction band. 
Thus neither C24B4 nor C24N4 were suitable candidates 
for reproducing the electronic structure of solid Cgo- 

The C28H4 molecule is found to be stable in the solid. 
Because the dangling bonds are passivated by hydrogen, 
apex-bonded C28H4 hyperdiamond (and other suggested 
structures 2 ^) are not favorable; the lower energy struc- 
ture is the face-bonded C28H4 hyperdiamond shown in 
Fig. 1(b). The latter is predicted to be a weakly bound 
solid with with a lattice constant of 16.3 A and a binding 
energy of approximately 0.2 eV per C28H4 molecule. The 
bonding between the six-membered rings is similar to 
that found for some orientations of Ceo molecules^ The 
structural and electronic properties of face-bonded C28H4 
hyperdiamond in both its pristine and doped forms are 
the subject for the remainder of this paper. 

Figure [21 compares the band structure and DOS be- 
tween solid Ceo and C28H4 around the band gap. The 
valence (conduction) band is formed from the three-fold 
degenerate LUMO (HOMO) of the C 28 H 4 molecule. The 
weakly broadened bands bear a striking similarity to 
those found in solid Cgo- Undoped solid C28H4 forms 
an insulator with 1-eV direct gap at T. Upon chemi- 
cal doping, the Fermi energy is expected to fall within a 
DOS peak with a value comparable to the alkali-doped 
Ceo materials. 

The dimensionless electron-phonon coupling parame- 
ter A = N(0)V cp depends on both the DOS at the Fermi 
energy iV(0) and the electron-phonon interaction poten- 
tial V cp which is proportional to the curvature of the 
fullerenei 3 i 7 i 8 That C28H4 is more curved than Ceo sug- 
gests that higher transition temperatures are obtainable. 
In general, V ep can be calculated as a double sum over 
the Fermi surface connecting states due to the phonon 
deformation potential. For molecular solids like Ceo 
and C28H4, the small dispersion of the electronic and 
phononic spectra implies that the molecular states and 
intramolecular phonons are an excellent approximation 
to those found in the solid. 2, 2j5, Therefore, we compute 
V ep for solid C28H4 by only considering the intramolecu- 
lar phonon coupling to the molecular states. 
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FIG. 2: Comparison of the band structure and DOS 
(states/eV/spin/cell) between solid fcc-Fm3 C60 and face- 
bonded C28H4 hyperdiamond. For C28H4 solid, the set of 
six bands above and below the gap are derived respectively 
from the three-fold degenerate LUMO and HOMO for each 
C28H4 in the two-molecule cell. The thin solid lines at zero 
show the centers of the direct gaps. 

The electron-phonon interaction potential is evalu- 
ated for the three-fold degenerate LUMO using the 
phonon frequencies ui v and eigenvectors e v for the iso- 
lated C28H4 molecule. Since our fullerene is composed of 
multiple species, it is convenient to include the mass of 
each species into the phonon eigenvectors normalization, 
Y^!i £ % u ' £ u'Mi — &w< , where N is the number of atoms in 
the fullerene. The electron-phonon interaction potential 
can be extracted from A 2 ^ and written in the form 

^ P = 4E^2 E I HE 4- v^K) I 2 , (i) 

" v v a.a' i 

where a and a' are the molecular states with degeneracy 
g for which the coupling is being evaluated. The ma- 
trix element in Eq. 0} is evaluated by a finite-difference 
approach where the deformation potential is calculated 
within the frozen-phonon scheme. Our calculated value 
of V ep is 181 meV, which includes contributions from 
most of the intramolecular phonons. It is observed 
in Raman-scattering experiments^^ and predicted by 
theorySiSZi that the A g phonons are screened out in alkali- 
doped Ceo- In the case of C28H4, the A g phonons con- 
tribute 12 meV to V cp . Even taking into account the 
screening of the A g phonons, the electron-phonon inter- 
action potential for C28H4 is over twice as large as that 
of C 60 (63 meV)J£ 

In Fig. 2 we see that analogous to Ceo, solid C28H4 is an 
insulator and can not superconduct unless doped. Ab ini- 
tio calculations of alkali-doped Ceo have demonstrated 
that intercalation of alkali atoms into the tetrahedral and 
octahedral sites result in the donation of the alkali va- 
lence electrons to the Cqq conduction band^SS* 2 ^^ This 
is reflected in the band structure by a lack of hybridiza- 
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tion of the alkali states with the conduction band, so that 
the band structure of the superconducting alkali-doped 
Ceo differs from that of pristine Cgo primarily by a rigid 
shift in the Fermi energy. This is viewed as the ideal 
doping case which we seek in the smaller fullerene solids. 

We have studied doping of solid C28H4 with the al- 
kali Na; other alkali atoms are expected to behave simi- 
larly. The intercalation of Na atoms into the solid C28H4 
was investigated in three different scenarios: (a) NaC2sH4 
intercalation into the tetrahedral site; (b) Na2C2sH4 
intercalation into the interstitial site between the six- 
membered rings on nearest neighbor C28H4 molecules; 
(c) Na@C28H4 endohedral doping. In each case, the Na- 
doped C28H4 structures were relaxed by conjugate gra- 
dient minimization of forces and stresses. We discuss the 
effects of doping as they pertain to the crystal structure, 
electronic properties and enthalpies of reaction^ 



(a) NaC 28 H 4 



(b)Na 2 C 2g H 4 



(c) Na@C 28 H 4 




FIG. 3: Band structure comparison among three doping cases 
(a) Intercalation of Na atoms into tetrahedral site; (b) Interca- 
lation of Na atoms into the interstitial sites between opposing 
six-membered rings on nearest neighbor C28H4; (c) Encapsu- 
lation of Na atom into the C28H4 cage. Dopant bands are 
depicted by dashed lines. The thin solid line at zero is the 
Fermi level. 

The intercalation of Na atoms into the tetrahedral site 
occurs with no significant expansion of the lattice [similar 
to that found, e.g., in K3C60 ( Ref. l3lh] . However, unlike 
the Ceo fullerides, additional bands appear just below the 
conduction band of C28H4 as depicted by the dashed lines 
in Fig.^a). A Mulliken population analysis reveals that 
the H atoms of the C28H4 gain additional charge, which is 
reasonable since four H atoms point into each tetrahedral 
site of solid C28H4 where a Na atom is located. Hence, in 
the NaC28H4 there is no doping of the conduction band; 
instead the H and Na states hybridize to form dopant 
states at the Fermi level. Furthermore, the calculated 
enthalpy of reaction indicates that the compound is not 
stable, so that this doping scenario would not be expected 
to occur naturally. 

Intercalation of Na atoms into the interstitial site be- 
tween opposing six-membered rings on nearest neighbor 



molecules results in a doping ratio twice as high as that 
of the previous case. Unlike NaC2sH4, there no longer 
are any dopant bands near the Fermi level. We deter- 
mine the enthalpy of reaction to be -0.9 eV per Na atom 
which is in line with those calculated for K3C60 3 ' 1 When 
the cell stress is relaxed, the lattice constant of Na2C2sH4 
increases from 16.3 A to 18.1 A resulting m a narrowing 
of the bands. In Fig.^b), we see the opening of a 0.2-eV 
gap at the Fermi level and so the solid transforms from 
a conductor into an insulator with a small gap upon cell 
relaxation. This is a direct result of Jahn- Teller distor- 
tion of the doped C28H4 molecules whose symmetry is 
reduced from Td to £>2d- The Jahn- Teller distortion of 
the isolated charged C28H4 molecule splits the 3-fold de- 
generate LUMO into a twofold degenerate and onefold 
degenerate state with a 0.2-eV gap. Hence, when the 
forces-and-stresses are relaxed in Na2C2sH4, a band gap 
similar in size to the Jahn- Teller gap opens up at the 
Fermi level. Thus in this case, doping by intercalation 
leads to insulating behavior. 

The encapsulation of Na atom into the C28H4 cage is 
the doping scenario which most closely parallels that of 
the alkali-doped Cgo solids. Figure [3{c) depicts a simple 
rigid shift of the Fermi level into the conduction band of 
C28H4. In contrast to Na2C2sH4, the lattice constant of 
Na@C28H4 does not increase. The Na atom is located at 
the center of the C28H4 and we find no Jahn- Teller dis- 
tortion. Although the isolated charged C28H4 molecule 
undergoes a Jahn- Teller distortion due to the degener- 
acy of the LUMO state, the broadening of the bands in 
Na@C28H4 is found to be sufficient to eliminate this ef- 
fect. Since it is known that endohedrally-doped C28 can 
be formed^ the relevant enthalpy of reaction for this 
system is the total energy of solid Na@C2sH4 relative 
to the isolated molecules; we find that this solid can be 
formed with an enthalpy of reaction similar to that of 
undoped C2sH4i2i 

Within the electron-phonon mediated theory of super- 
conductivity, A = N(Q)V CV plays a crucial role in deter- 
mining T c . The superconducting transition temperature 
for Na@C28H4 can be estimated using McMillan's solu- 
tion of the Eliashberg equations, 32-33 



1.2 



exp 



1.04(1 + A) 



A- M *(1 + 0.62A) 



(2) 



where oj\ n is a typical phonon frequency and /i* is the 
Coulomb pseudopotential which describes the effective 
electron-electron repulsion. Typical values of uj\ n ~ 
10 3 K for C28H4 and Ceo- We may expect that /x* 
for C28H4 should not differ considerably from that of 
Ceo since the subbands and phonon energies are sim- 
ilar. For alkali-doped Ceo solids, experimental results 
lead to n* « 0.22.24 The endohedral doping scenario 
for C28H4 gives that iV(0) « 5 states/eV/spin/molecule 
which is half of that found in the canonical alkali-doped 
Ceo structure, i.e., K 3 Cqq^ Accounting for screening 
of the A g phonons in the electron-phonon interaction, 



a I4i 



169 meV gives an enhancement factor of 
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A(Na@C 28 H 4 ) w 1.5A(K 3 C 60 )— These approximations 
are in line with those found in the literatureSiiii& and is 
sufficient for qualitative comparison. Using this value of 
H* and wi n , one finds that A(K3C@o) = 0.84 in order to ex- 
plain the experimentally observed ^(K^Ceo) = 19.3 K. 35 
The enhancement factor for A that we calculate predicts a 
T c (Na@C2 8 H 4 ) « 3T C (K 3 C 60 ) « 58 K, higher than that 
found in the highest temperature alkali-doped fulleride 
[40 K in Cs 3 C 60 (Refill]. 

In summary, ab initio pseudopotential density func- 
tional calculation were performed on C2 8 -derived molec- 
ular solids. We find a C28H4 solid which binds weakly 
and exhibits many of the band structure features as Cgo- 
Several doping scenarios are investigated and we find that 
Na@C28H4 to be the most promising for superconductiv- 
ity. The calculated electron-phonon interaction potential 
and DOS leads to a T c enhancement of three times that 



found in K^Ceo- Since endohedral doping is expected to 
be insensitive to the type of dopant, this suggests promis- 
ing possibilities with other atoms including ones of higher 
valence that can control the doping level. Guo et alm& 
have succeeded in encapsulating group-IVB atoms in C28, 
which suggests that other small fullerenes such as C28H4 
can also be endohedrally doped. 
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